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AERONAUTICAL SYMBOLS 
1. FUNDAMENTAL AND DERIVED UNITS 



y 


Symbol 


Metric 


English 


Unit 


Symbol 


Unit 


Symbol 


Length 


I 

t 

F 


meter 


m 
s 

kg 


foot (or mile) 


ft. (or mi.) 
sec. (or hr.) 
lb. 


Time 


second , 


second (or hour) 


Force 


weight of one kilogram 


weight of one pound 


Power 


P 


kg/m/s 




horsepower 


hp », 
m. p. h. 

f. p. s. 


Speed 


fkm/h 


k. p. h. 
m. p. s. 


mi./hr 






W/s 


ft./sec « 



W, Weight = mg 

g, Standard acceleration of gravity = 9.80665 
m/s 2 = 32.1740 ft./sec. 2 

m, Mass = — 



p, Density (mass per unit volume). 
Standard density of dry air, 0.12497 (kg-m"" 4 

s 2 ) at 15° C. and 760 mm = 0.002378 

(lb.-ft.- 4 sec. 2 ). 
Specific weight of "standard" air, 1.2255 

kg/m 3 = 0.07651 lb./ft. 3 . 



GENERAL SYMBOLS, ETC. 

mk 2 y Moment of inertia (indicate axis of the 
radius of gyration k y by proper sub- 
script). 
Area. 

Wing area, etc. 
Gap. 
Span. 
Chord. 



S, 

G, 
b, 

c, 
b 2 
S' 

Mi 



Aspect ratio. 
Coefficient of viscosity. 



V, 
2. 
L, 
D, 



True air speed. 



3. AERODYNAMICAL SYMBOLS 

Q, Resultant moment. 
Q 



Dynamic (or impact) pressure =^pV 2 . yj 

Lift, absolute coefficient C^^U 

q£> 



Drag, absolute coefficient C D - 



D 

"qS 



Do 



Di, Induced drag, absolute coefficient C Di 
force, absolute 



C, 



R, 

I'M) 

it, 



Profile drag, absolute coefficient @ D o = ^j 

Dj 
qS 

Parasite drag, absolute coefficient C Dp =^ 
Cross-wind 

Cc ~qS 
Resultant force. 

Angle of setting of wings (relative 

thrust line). 
Angle of stabilizer setting (relative 

thrust line}. 



qS 

coefficient 



to 



to 



Resultant angular velocity. 

t Reynolds Number, where I is a linear 

dimension, 
e. g., for a model airfoil 3 in. chord, 100 
mi./hr. normal pressure, at 15° C, the 
corresponding number is 234,000; 
or for a model of 10 cm chord 40 m/s, 
the corresponding number is 274,000. 
C py Center of pressure coefficient (ratio of 
distance of c. p. from leading edge to 
chord length). 
Angle of attack. 
Angle of down wash. 
Angle of attack, infinite aspect ratio. 
Angle of attack, induced. 
Angle of attack, absolute. 

(Measured from zero lift position.) 
Flight path angle. 
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A COMPARISON BETWEEN THE THEORETICAL AND MEASURED LONGITUDINAL 
STABILITY CHARACTERISTICS OF AN AIRPLANE 

By Hartley A. Soul£ and John B. Wheatley 



SUMMARY 

This paper covers an investigation of the application 
of the theory of dynamic longitudinal stability, based on 
the assumption of small oscillations, to oscillations an 
airplane is likely to undergo inflight. The investigation 
was conducted with a small parasol monoplane for the 
fixed-stick condition. The period and damping of longi- 
tudinal oscillations were determined by direct measure- 
ments of oscillations inflight and also by calculation in 
which the factors that enter into the theoretical stability 
equation were determined in flight. A comparison of 
the above-mentioned characteristics obtained by these two 
methods indicates that the theory is applicable to the 
conditions encountered inflight. 

The investigation was extended to determine the feasi- 
bility of calculating the stability characteristics from basic 
data of the type that would be available to a practicing 
designer. The results of this phase of the work show that 
for the power-off condition the agreement between the 
actual and predicted stability characteristics was reason- 
ably satisfactory, except perhaps when the predicted sta- 
bility is close to neutral. For the power-on condition the 
stability can not be predicted owing to the present lack of 
information concerning slipstream effects. Further prog- 
ress in a solution of the problem of longitudinal stability 
depends largely on increasing the knowledge of these 
effects. 

INTRODUCTION 

With each new airplane design, the designer is con- 
fronted with the problem of attaining some, although 
not a definitely established, degree of longitudinal 
dynamic stability. The theory utilized in this prob- 
lem has been evolved and published in various stand- 
ard works such as those of Bryan, Bairstow, Cowley 
and Levy, and Glauert. This theory involves the 
basic assumption of small deviations from the steady 
state and the application of numerous data relating 
to the airplane. The theory is not widely applied in 
design, probably because an acceptable (although not 
necessarily the most desirable) degree of dynamic sta- 
bility is often attained by the application of simple 
rules developed from experience, and also because the 
theory is fairly complex and its validity has not been 
clearly demonstrated. New designs, however, are fre- 



quently found to possess decidedly objectionable 
stability characteristics. Furthermore, it is quite 
probable that the stability characteristics of many 
airplanes could be considerably improved. Owing to 
the lack of knowledge concerning the factors that tend 
to produce dynamic stability and the relative impor- 
tance of these factors, the nature of the changes nec- 
essary for improvement is not generally understood. 
Thus there has arisen the need for a comprehensive 
study of the problem to ascertain the validity of the 
theory, to determine in what respects, if any, it requires 
modification, and to establish thereby the procedure 
by which the designer can predict the stability char- 
acteristics of an airplane from basic data. An addi- 
tional phase of the problem arises from the fact that 
the most desirable degree of stability is not definitely 
established. The solution to this phase probably 
should be obtained by determining the degree of 
dynamic longitudinal stability possessed by airplanes 
that have very good flying characteristics. 

For the reasons mentioned above, an investigation 
of the dynamic longitudinal stability of an airplane in 
flight with and without power was undertaken. The 
results of the investigation, as reported herein, are 
divided into two parts. The first deals with checking 
the validity of the basic assumption of small devia- 
tions from the steady state against the actual condi- 
tions encountered in flight. The second deals with 
the calculation of stability from basic data and the 
comparison between calculated and test results. The 
flight tests required by the investigation were made 
with a small parasol monoplane. No attempt was 
made in the present case to establish the most satis- 
factory degree of dynamic longitudinal stability. 

THEORY OF STABILITY 

GENERAL PRINCIPLES 

As previously mentioned, the classical theory of the 
stability of small oscillations, its application to the air- 
plane, and the development of the equation of dynamic 
longitudinal stability are completely covered in the 
various works on the subject. The standard form of 
the stability biquadratic, however, is given in dimen- 
sional units and is very difficult to handle. Glauert 
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I lis developed a nondimensional form for the equation 
(reference 1) which is more convenient in application 
and which greatly simplifies the analysis of the longi- 
tudinal stability. This nondimensional form has been 
used in the present case. 

For convenience of reference, the transformation 
from the dimensional to the nondimensional form of 
the stability equation is given in the appendix using 
the absolute system of units. As explained therein, 
only the long-period, lightly damped oscillation is con- 
sidered. The nondimensional expression for the period 
of this oscillation is 



IE ( DC— BE Y 



and that for the damping coefficient is 



fl ~ 2\G &) 



where 



B= -m q -x u -z w 

0= z w m q + ZwX u + m q x u — z u x w - \xm w 

D = ^nm u C L + fjtm w x u + ^ C L tan 8 0 im w 

+ m q (z u x w — x u z w ) — x w \im u 

and 

E= \ v>Cl tan O 0 (x w m„ - x u m w ) + ^ p-C L {m w z u - m u z w ) 

If the damping coefficient is negative, the motion is 
stable; that is, the oscillations are damped; and if it is 
positive the motion is unstable. 

The terms x u , x m z U) z m m u , m m and m q are the non- 
dimensional derivatives (or derivative coefficients), 
and /x is a parameter. As noted in the appendix the 
length I on which m u , m w , m Q) and p. are based can be 
chosen arbitrarily. The length used throughout this 
paper is the distance from the center of gravity to the 
rudder post. The terms C L and 8 0 refer to the initial 
condition. 

The above expressions for the period and damping 
can be readily converted to the dimensional form for 
comparison with experimental data by introducing the 
pVS 



factor 



and 



m 



T ~ Tl P VS 



IK 



m 

where T and f are the period and the damping factor, 
respectively, in dimensional form. 

THE RESISTANCE DERIVATIVES IN TERMS OF FUNDAMENTAL AIR- 
PLANE CHARACTERISTICS 

The terms x u , x w , C L , 0 O , z u , z w , m«, m„, m„ and the 
parameter n must be evaluated for an analysis of the 
dynamic longitudinal stability. The manner in which 



most of the derivatives can be expressed in terms of the 
fundamental airplane characteristics will now be shown. 
Because of the marked change produced by the pro- 
peller in the stability derivatives, the derivatives will 
be first discussed for the power-off condition, then for 
the power-on condition. 

Power-off. — -A basic assumption for this case is that 
with the propeller idling the V/nD ratio for the propeller 
remains constant for small changes in V. The assumed 
initial conditions is a steady glide with the propeller 
operating at the V/nD of zero effective thrust. 

In the appendix it is shown that for the power-off 
condition x w =—C D ; similarly, it can be shown that 

z u = —C L and m u = ^j^- 

These expressions for x u and z u apply when C D and 
C L are taken as the over-all drag and lift coefficients of 
the airplane. As the airplane is in the steady gliding 
condition C m is zero and consequently m u = 0. 

It is noted at this point for later reference, that in 
actual gliding flight with throttled engine the pro- 
peller, in general, develops a small effective thrust, 
T 7 , which may have an appreciable effect on the term 
Xw In this case 



X- 



T-D 



C D %pV<S 



and for constant V/nD 



Xu 



where C D ' is the effective drag coefficient of the air- 
plane-propeller combination. 

The velocity w introduces a change in the magnitude 
and direction of the resultant velocity. The change in 
magnitude, however, is a second order effect and can be 
neglected. The change of angle causes a change in 
the flight-path angle. Then 



dw Aa Ann 



dZdtf 
d0 dw 



tan0 = 0 = 7 and^=y 

for a small angular displacement 6 

X= - C D ^S cos 0+ C l *y-S sin 6 

as 6 is small, cos 6=1 and sin 8=6 

Also, it may be assumed that 6 equals the change in a. 

Then 

Xw ~ ' 2 da +/20l/ 
By a similar procedure it is found thai 
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The term m w is given by the expression 



, C dC m 

2 l V da 



The coefficient m Q of the moment derivative M q is a 
basic characteristic of the airplane comparable with 
C L , C d , and C m . The expression for this coefficient is 
simply 

m q = MJpSWr) 

The term 11 M q " can be determined from calculations 
by the method described later in the paper, or by direct 
measurement with a model or full-sized airplane. 

The foregoing discussion shows that a complete 
theoretical analysis of the dynamic longitudinal 
stability for the power-off condition can be readily 
made from basic data. The basic data required for a 
complete determination of the power-off stability 
characteristics can be summarized as follows: 
Wing area. 
Mass of the airplane. 

A characteristic length, I (here taken as the 
distance from the center of gravity to the 
rudder post). 

Moment of inertia about the Y axis. 

Lift, drag, pitching moment, and rotary deriva- 
tive coefficient as functions of angle of attack. 
Power-on. — As already pointed out, the propeller, 
because of its thrust and slipstream, causes consider- 
able change in the stability derivatives. In the first 
place the propeller, when operating at a positive 
thrust, must have a distinct series of stability deriva- 
tives of its own which are additive to those of the 
airplane. (Reference 2.) It is very probable that the 
majority of these, particularly those depending on 
angular velocity, are negligible. The variation of 
longitudinal force or thrust with velocity, however, 
must be considered. Besides the derivatives of 
the propeller, there is a large change produced in the 
stability derivatives of the airplane because of rela- 
tively high velocity of the slipstream flowing over the 
central portion of the wings and over the fuselage and 
tail. The increase of drag due to this effect need not 
be considered when the propeller characteristics are 
based on effective thrust. 

Our present knowledge of the character of the flow 
behind a propeller does not permit a satisfactory pre- 
diction of the effect of the slipstream on the stability, 
although the experiments discussed later show the 
magnitude of the effects on some of the derivatives for 
this particular airplane. The following discussion of 
the power-on derivatives will be carried as far as our 
present knowledge of the propeller effects allows. 

A basic assumption in this case is that the propeller 
torque remain constant. The assumed initial condi- 
tion is steady level flight. 



The thrust now varies with velocity so that we have 



1 AT 
pSVdV 



-C D 



As T is here taken as the effective thrust, C D is the 
same as for the power-off condition. An expression for 
A T 

in terms of generalized propeller characteristics has 

been derived by Glauert (reference 2) for the assumption 
of constant torque. This expression is 



AT 



where 



27 V (2a-3r + r 3 ) 
" V (1-r 2 ) (2a + r 3 ) 

a is a constant = 1.325 



and r is the ratio of V/nD of the initial condition 
to V/nD for zero thrust. 

Consideration of Glauert 's expression shows that, 
d T 

although the term -^y varies considerably with changes 

in r over the normal speed range, the term is always 
small relative to C D . Consequently, it is permissible to 
evaluate the term on the basis of the average value of r. 
Then we can write 



dV V /v 



1 dT 
P SVdV 



since the thrust equals the drag when the angle of 
attack of the thrust line is neglected. Then 

The form of the expression for z u is the same as 
for the power-off but, owing to the previously men- 
tioned slip stream effect, the lift coefficient at a given 
angle of attack is in reality not the same as in the 
power-off condition. Actually 

z u =-C L '=-(C L + AC L ) 

where AC L represents the effects of the slipstream and 
the thrust component on the lift coefficient. The term 
AO L can not be satisfactorily evaluated. The experi- 
ments discussed later show that at least in this particu- 
lar case this term is not negligible. 

The principal difficulty in the evaluation of the 
power-on stability derivatives pertains to the moment 
derivatives. Both the thrust and slipstream effects 
enter into these terms. The thrust moment is intro- 
duced because in the general case the line of action of 
the propeller does not pass through the center of 
gravity. The action of the slipstream is more complex. 
The change of lift coefficient of the wings is accom- 
panied by a change of downwash angle which in turn 
produces a change of angle of attack at the tail. In 
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addition to the change of angle of flow at the tail, the 
velocity there is increased. As the variation of the 
angle and velocity at the tail with a change in the 
airplane's velocity is not known, it is impossible to 
evaluate m u . That m u is not zero for the power-on 
condition is demonstrated by the change of balance 
that takes place when power is applied and the V/nD 
of the propeller changed from that for zero thrust to 
some other value. The factors that complicate the 
solution of m u also enter into m w . The moment deriva- 
tive m q is less complex and varies from the power-off 
derivative essentially because of the relatively greater 
velocity over the tail for the power-on condition. 

Attempts were made to evaluate the derivatives 
m M) m w , and m q for the power-on condition on the 
basis of assumed flow over the tail. Although reason- 
able agreement was obtained for the damping and period 
calculated from computed values of these derivatives, 



from the type of data obtainable before the airplane 
is constructed. Complete measurements were made 
for the power-off condition. For the power-on condi- 
tion only the derivatives depending on C L and G D 
were determined as it is impossible to measure the 
moment derivatives for this condition. 

THE AIRPLANE 

The airplane used in the investigation was a Doyle 
0-2, a small parasol monoplane. (Figs. 1 and 2.) The 
dimensions required in the study of longitudinal 
stability are given in the following table: 

Weight 1,290-1,340 lb. 

Moment of inertia about Y axis (B) 673 slug ft. 2 

Wing dimensions: 

Span 30.0 ft. 

Chord 5.50 ft. 

Area 159.5 sq. ft. 

Airfoil section Clark Y. 




Figure 1.— Front view of Doyle 0-2 airpl 



so many arbitrary assumptions had to be made in 
their evaluation that it is not considered worth while 
to include either the methods used to compute the 
derivatives or the results of the calculations made 
with them. 

EXPERIMENTAL STABILITY 

Flight tests were made with the airplane utilized 
in this investigation to determine the stability 
derivatives previously discussed and also to determine 
the period and damping of oscillations in flight by 
direct measurements. These data served two purposes. 
By comparing the period and damping found by 
direct measurements with the period and damping 
calculated from derivatives obtained from flight 
data, it was possible to test the applicability of the 
theory. Also, by comparing the derivatives and sta- 
bility characteristics obtained from flight data with 
those calculated from basic data for the airplane 
elements, as explained later, it was possible to test 
the precision with which stability can be predicted 



Ccnter-of-gravity position: 

Below chord 23.4 in., 36.0 per cent 

chord . 

Rear of L. E 22.4 in., 34.0 per cent 

chord. 

Above thrust line (h) 12.0 in. 

Horizontal tail surface: 

Span 9.0 ft. 

Area 18.7 sq. ft. 

Airfoil section Flat plate. 

Distance from c. g. to rudder post (I) 11.8 ft. 

Engine LeBlond, 60 hp 

Propeller pitch-diameter ratio 0.60 

Wing loading 8.09-8.40 lb./sq. ft, 

Power loading 21.5-22.3 lb./hp 

INSTRUMENTS 

The instruments used during the flights were : 
Air-speed recorder. 
Recording inclinometer. 
Timer. 

Control-posi ti( m recorder. 
Angular-velocity record e r . 
Angle-of-attack recorder. 
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All the instruments are photographic recording. 
Detailed descriptions of the air-speed recorder, in- 
clinometer, timer, control-position recorder, and angu- 
lar-velocity recorder are given in references 3, 4, 5, 6, 
and 7, respectively. The angle-of-attack recorder is a 
recent development. It consists of a light wind vane, 
free to rotate about an axis parallel to the Y axis of the 



angle was obtained from the angle of attack and incli- 
nation of the X axis. The lift and effective drag 
coefficients were deduced from the wing area, dynamic 
pressure, and components of weight parallel and per- 
pendicular to the glide path, the effective drag coeffi- 
cient being the drag coefficient of the airplane-propeller 
combination. In addition to the above-mentioned 




Figure 2.— Side view of Doyle 0-2 airplane 



airplane, mounted about a chord length ahead of the 
wing at the outboard strut fitting. The motion of the 
vane is transmitted through the supporting tube by a 
fishline to the recording mechanism. The vane and 
recording mechanism are so placed in relation to one 
another that when the vane is in the desired position 
the recording mechanism may be embedded in the 
wing so as not to disturb the flow. A schematic 
sketch of the angle-of-attack recorder is shown in 
Figure 3. 



FLIGHT PROCEDURE 



Preliminary flights. — Preliminary to the main flight 
tests it was necessary to calibrate the air-speed system 
and the angle-of-attack recorder in flight. For this 
purpose a trailing Pitot-static tube, a second air-speed 
recorder, and an indicating statoscope were added. 
The air-speed system was calibrated against the trail- 
ing Pitot-static tube suspended 70 feet below the wing. 
The angle-of-attack recorder was calibrated against the 
inclinometer in the airplane during steady level flight . 
The calibrations were made over the complete speed 
and angle-of-attack range. 

Determination of C L and C D * — Lift and drag co- 
efficients for the power-off condition were determined 
in glides with the throttle closed. The dynamic 
pressure, angle of attack, and inclination of the X 
axis were recorded directly in glides. The gliding 



items, the average altitude, engine speed, and air 
temperature were noted during each glide. 

Lift coefficients for the power-on condition were 
obtained by recording dynamic pressure and angle of 
attack in a series of level-flight runs. The true drag 
coefficients for this condition were deduced from the 
results of the glide tests by correcting the effective < I rag 
coefficients for the calculated drag of the idling pro- 



Film 
drum" 




! o j Recording instrument 
Wing contour 

'■ b , Mirror Support tube 



Cord for transmitting--' 
motion of wind vane to 
recording instrument 




Wind vane 



Figure 3.— Schematic sketch of the angle-of-attack recorder 

peller. The thrust coefficients required in these 
calculations of propeller drag were determined from 
the actual V/nD of the propeller in each glide and a 
curve of calculated propeller characteristics. 

Determination of C m , — The static pitching moment 
was obtained from the results of glide tests in which a 
known pitching moment was applied. The pitching 
moment was applied by a sliding weight mounted so 
that its position could be controlled from the pilot's 
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cockpit, A 1 00-pound lead weight was mounted on a 
steel tube in the front cockpit and was fitted with a [ 
control and locking device so that it could be locked at 
the center of gravity and at various known distances 
from the center of gravity. (Figs. 4 and 5.) 

The airplane was first glided with the weight at the 
center of gravity at a given air speed and with the 
elevator held stationary. Records were made of the 
air speed, angle of attack, and elevator position. The 
elevator position was recorded to make certain that 
the elevator angle remained constant during the runs. 
Subsequent records were made with the weight moved 
to a new position and with the original elevator angle. 




Figure 4.— The installation of the sliding weight in the front cockpit 



Determination of m q . — The damping coefficient due i < » 
rotation for the power-off condition was determined by 
studying the motion of the airplane during pitching 
oscillations in glides. The pitching oscillations were 
set up from an initially steady state by a fore-and-afi 
movement of the stick after which the stick was re- 
turned to its original position. Continuous records of 
the angular velocity in pitch, air speed, angle of 
attack, and elevator angle were obtained during a por- 
tion of the steady glide and an ensuing period of several 
seconds. Angular acceleration was deduced from re- 
corded angular velocity by graphical differentiation. 
The resultant pitching moment acting on the airplane 




Figure 5. — The pilot's cockpit showing crank for moving the sliding weight 



This procedure was repeated until the entire normal at any given instant was calculated from the accelera- 

range of angle of attack with one elevator position was I tion and the moment of inertia of the airplane about 

covered. Difficulty experienced in maintaining a con- the Y axis. The pitching moment caused by the 

stant elevator angle was obviated by providing a stop rotation was found by deducting the static component 

on the control column. from the resultant and was then reduced to coefficient 

The static pitching moments were calculated from form. The static component of the resultant pitching 

the weight and position of the sliding mass reduced to moment for a given instantaneous angle of attack was 

coefficient form and plotted against angle of attack. ! found by referring to the angle of attack in the initial 

The moment-coefficient curve was obtained for only I steady state, the curve of static pitching-moment 

one elevator setting. In the application of these j coefficient against angle of attack previously estab- 

results to the calculation of stability the assumption lished, and the recorded air speed for the same instant, 

was made that a change of elevator angle produces a As the static pitching moment is zero in a steady glide, 

constant difference in the ordinates over the entire the reference to the angle of attack of the initial steady 

extent of the curve. condition provided a means of determining the shift 
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in ordinates of the curve of static moment coefficient 
against angle of attack corresponding to any given 
elevator position. The static pitching-moment coef- 
ficient for any given elevator position and angle of 
attack was thereby readily established. Values of m q 
for the power-off condition were calculated in the 
above-described manner for various angles of attack 
throughout the range of normal flight. 

Determination of period and damping coefficients. — 
The period and the damping coefficient were deter- 
mined by direct measurements of the oscillation char- 
acteristics in flight at the same altitude as that at 
which the previously mentioned tests were made. 
These direct measurements were made both for the 
power-off and power-on conditions. As u, v)> and 6 
are interdependent variables, the periods of their va- 



whore 

f, damping factor. 
T, period, in seconds. 

Vi, V 2 , V z , true air speed in feet per second at two 
successive maximums and the intervening mini- 
mum. 

PRECISION OF MEASUREMENTS 

Frequent check calibrations were made of the instru- 
ments used in the flight tests. Appreciable errors 
caused by faulty operation or changed calibration of 
the instruments were thereby eliminated. The effect 
of lag in the angular- velocity recorder was eliminated 
by a correction determined in laboratory tests. Prob- 
ably the most serious source of error in the results was 




Figure G.— Typical records of air speed during oscillations 



riatikms with time are necessarily the same, although 
tike variations may not be in phase with each other. 
( Consequently, the period and damping could be deter- 
mined by studying the behavior of only one variable. 
Air speed was chosen as the one most convenient for 
study. The tests were made by setting up oscillations 
in the manner previously described and obtaining 
records of air speed for at least one complete oscilla- 
tion. Sample records are shown in Figure 6. The 
records designated Runs 1 and 2 illustrate an unstable 
condition experienced with power-on at initial angles 
of attack of 16.5° and 15.0°, respectively; and those 
designated Runs 3 and 4 illustrate a stable condition 
with power-off at an initial angle of attack of 4.8°. 

The period of the oscillation was found by the time 
interval between two consecutive points of maximum 
velocity. The damping coefficient is given with suffi- 
cient accuracy by the approximate relation 



introduced by the limitations of accuracy in determin- 
ing angular accelerations by graphical differentiation 
of the angular- velocity records. Wide record lines 
caused by instrument vibration were a contributing 
cause to inaccuracies in evaluating records. The 
effect of accidental errors was* considerably reduced, 
however, by obtaining two or three repeat runs for 
each test condition and expressing the results as faired 
curves. Following is a list showing the estimated 
limits of accuracy for the various items measured and 
derived from flight tests: 

Experimental quantity Precision. 

C L (except at C Lm ax) ± 2 per cent to ± 5 per cent. 

C L maz ±10 per cent. 

C D ±2 per cent to ± 8 per cent. 

C m ±5 per cent. 

m Q ±20 per cent. 

__ __1 . xl f±10 per cent. 

Period../ e*Peri*"ntal... "{±5 per cent. 
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CALCULATED STABILITY 

In most cases it is desirable to complete an analysis 
of the stability of a proposed airplane design before 
construction. For such an analysis there are at hand 
standard basic data regarding the aerodynamic char- 
acteristics of the airplane elements from which the 
aerodynamic characteristics of the complete airplane 
can be calculated. In an advanced stage of the design 
there may also be available results of wind-tunnel 
tests on a model of the airplane from which the flight 
chiiracteristics of the airplane can he more accurately 
predicted than from data for the airplane elements. 
In the present case an analysis of the dynamic longi- 
tudinal stability of the airplane used in the flight tests 
has been made by utilizing only basic data for the 
airplane elements; i. e., wing airfoil characteristics, 
parasite drag of nonlifting elements, and tail-plane 
characteristics neglecting interference effects and 
assuming a tail efficiency of 1. 

It will be assumed at the outset that the design has 
been carried to such a point that certain principal fea- 
tures of the airplane have been tentatively adopted. 
These features include the weight, the wing profile, 
area, dimensions, and arrangen cut with respect to the 
center of gravity; the horizontal tail-plane profile, area, 
dimensions, and position relative to the center of 
gravity and to the wings; and the over-all length and 
height, including the landing gear. In order to study 
the longitudinal stablility it is then necessary to find in 
addition the moment of inertia about the Y axis and 
the lift, drag, static pitching-moment and rotary 
pitching-moment coefficients as a function of angle of 
attack. 

The moments of inertia of the airplane used in 
these tests had previously been measured (reference 8); 
consequently, the moment of inertia about the Y axis 
was available for the calculations. In general, this 
moment of inertia will not be known and must be 
estimated. It may be determined from the summation 
of small elements. This method, however, besides 
requiring long and tedious calculations, is not likely 
to prove very accurate. Probably the most convenient 
and accurate method is that given in reference 8 in 
which there are given the results of accurate measure- 
ments of the moments of inertia of several airplanes 
representative of different wing arrangements. From 
these results there have been derived nondimensional 
coefficients based on the weight and over-all dimen- 
sions. The suggested procedure is to choose the co- 
efficients of the airplane most nearly similar to the 
projected design and to obtain the moments of inertia 
of the projected design by multiplying the nondimen- 
sional coefficients by the appropriate factors. 

Curves of C L and C D for the complete airplane were 
obtained by consideration of the airfoil character- 
istics as determined in tests at large scale in the vari- 
able-density w T ind tunnel, and of the probable effect 



of the other airplane elements. The lift coefficients 
for the airfoil, after correction to the aspect ratio of 
the actual wing, were assumed to apply to the complete 
airplane. To the drag coefficients of the airfoil were 
added an estimated drag coefficient representing the 
effect of the fuselage and engine, landing gear, and 
exposed struts. The parasite-drag coefficient was 
assumed to remain constant at all angles of attack. 

The static pitching-moment coefficient was obtained 
by first constructing a vector diagram of the wing 
resultant force at a constant air speed. On the dia- 
gram the center of gravity of the airplane was spotted 
in its correct position in relation to the wing chord. 
The moment of the wing about the center of gravity 
was then computed with the data obtained from the 
diagram. As the moment of the parasite drag about 
the center of gravity was assumed to be zero, the only 
other moment was derived from the tail The down- 
wash angle at the tail was computed from the equations 
developed b}^ Diehl (reference 9) employing the con- 
stant recommended by Reid (reference 10) for general 
use, 

e = 60//? Cr + l)-°- 38 (y+l)-°- 2Z C L 

where x and y are the horizontal and vertical distances 
of the (ail plane behind and below the wings in units of 
chord length. The angle of attack of the tail was 
corrected for the downwash from the wing. Next the 
slope of the lift-coefficient curve for the tail plane was 
estimated by consideration of the airfoil section and 
aspect ratio. The tail moment about the center of 
gravity was then calculated for the same velocity as 
for the wing, from the slope of the tail lift-coefficient 
curve, angle of attack of the tail, the tail area, and the 
moment arm. This moment w r as added algebraically 
to the moment of the wing and the resultant moment 
reduced to coefficient form. 

The coefficient m Q of the rotary derivative was cal- 
culated directly from the tail characteristics. The 
assumption was made that the damping of rotation of 
the tail is not materially augmented by the damping of 
the other parts of the airplane and that the total damp- 
ing can be taken as equal to that of the tail. The 
equation used for the calculation of m q may be derived 
as follows: Let the change in angle of attack of the tail 
caused by the angular velocity q be &a T 

Aa r = w T /V T 

now 

w T = Iq 

and 

v T ={v*+(i q y}» 

Let M T be the change in moment of the tail caused by 

the angular velocity q 

then 

dC 

M T = l*q-^ S T %{V 2 +(lqY)' A 
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as lq is small compared to V the term (lq) 2 may be 
neglected 

dM 7 



and 



' "2"^~d^r 



M Q = y S T 1 dC LT 
a ~pVSl 2 rj h S v da T 



REvSULTS AND DISCUSSION 

The results of the investigation are given in Figures 
7 to 16, inclusive. The lift and drag curves of the air- 
plane are given on Figure 7, showing the power-off 
curves obtained both from flight tests and from cal- 
culations. The experimental power-on curves of lift 
and drag are given for comparison. In Figure 8 the 
experimental and calculated curves of pitching-moment 
coefficients are plotted. Figures 9, 10, 11, 12, and 13 
show the stability derivatives computed from the data 
of the two previous figures, and Figure 14 gives the 
experimental and calculated curves of the stability 
derivative m Q . In Figure 15 are plotted curves of 
period and damping obtained from the experimental 
and calculated stability derivatives and points repre- 
senting the actual measured period and damping. 
Faired curves are given in Figure 16 to show the differ- 
ence in the measured stability characteristics power- 
oil' and power-on. It should be noted that all results 
are given for an air density of 0.00217 slug per cubic 
foot which corresponds to an altitude of 3,000 feet in a 
standard atmosphere. 

For the determination of the validity of the assump- 
tion of small oscillations in dealing with the stability 
it is obviously necessary that the disturbances during 
which the actual period and damping were measured 
should have been large enough to represent actual 
average conditions that the airplane is expected to 
meet in normal flight. An inspection of the air-speed 
records obtained during the oscillations showed varia- 
tions in magnitude as great as 30 miles per hour during 
the course of several of the runs and average deviations 
in excess of 20 miles per hour. The amplitudes of 
these oscillations are thought to be great enough to be 
regarded as representative of conditions following a 
disturbance produced by the roughness of the air. 

The nature of the agreement between the theory 
and experiment is shown on Figure 15 where the curves 
of period and damping coefficients as calculated from 
the experimental derivatives are shown with the 
points representing the measured values of the same 
items obtained from the oscillation tests. It can be 
seen that although none of the points actually fall 
on the theoretical curve of damping coefficient they 
are dispersed about equally on either side of the curve. 
The dispersion of the points is probably caused to a 
large extent by unsteady air conditions affecting the 
airplane during the time the records were made. 
The agreement obtained indicates that the theory 



is satisfactory for use in analysis of the longitudinal 
stability characteristics of airplanes. 

For the prediction of stability characteristics it is 
essential that not only can the stability characteristics 
be calculated from the derivatives but also that the 
derivatives can be calculated with the necessary degree 
of accuracy. The degree of accuracy with which the 
stability derivatives can be calculated is shown by a 
comparison of the calculated and experimental curves 
for the power-off condition in Figures 7 to 14. The 
curves of lift coefficient and associated derivatives 
show fair agreement except at the stall. The calcu- 
lated drag curve differs from that obtained in gliding 
flight by an amount that is approximately equal to the 
drag attributed to the idling propeller and conse- 
quently is in fair agreement with the curve for the 
power-on condition. This discrepancy, which is 
brought about by the neglect of the drag of the idling 
propeller, is reflected in x u and x w where the disagree- 
ment is proportional to that for C D . The greatest dis- 
crepancy between the calculated and experimental 
derivatives is in C m and m w and is probably caused by 
a nonelliptical lift distribution and interference effects 
at high angles, with the result that the calculated 
downwash angles are in error. With m g the agreement 
is fair. 

The results of the computations of period and damp- 
ing from the calculated derivatives are of interest 
notwithstanding the relatively poor agreement in 
ra w which seemingly has no serious effect on the final 
result. The curves of the calculated period and 
damping coefficient show the same general trend as 
those computed from the experimental derivatives, 
though the damping curve is offset an appreciable 
amount. The calculated values for the damping are 
smaller than the true values for this case. The com- 
parison indicates that predicted results would often 
be useful where more exact data are lacking. Without 
the assurance of better accuracy than that attained in 
the present case, however, it seems evident that 
predicted results showing close to neutral stability 
should not be interpreted too literally. In this con- 
nection it should be noted that much better accuracy 
can probably be attained if model test data are avail- 
able for the analysis. 

Aside from the principal results of the report there 
are other points of interest that have been brought out. 
The most important of these is the marked difference 
in the stability characteristics between the power-off 
and power-on conditions. (Fig. 16.) The stability 
for the power-on condition is generally less than with 
power-off. A portion of this difference can be attrib- 
uted to the difference in flight path. Supplementary 
calculations made show, however, that the difference 
in flight path accounts for very little of the difference 
and that the larger part must be attributed to slip- 
stream effects. It will be further noticed that the 
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curvature of the power-on damping coefficient curve 
is reversed at high angles of attack and with a very 
slight change of angle the damping coefficient becomes 




O 4 8 12 16 20 24 

oi, angle of attack, degrees 

Figure 7. -Lift and drag coefficients of Doyle 0-2 airplane 
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Figure 8.— Static pitching-moment coefTicient of Doyle 0-2 airplane 

positive, i. e., the motion becomes unstable. This 
change of curvature can not be explained by making 
allowance for the effect of the thrust and must be 
dependent on the slipstream effect. Another demon- 
stration of the slipstream effect is shown on Figure 



7 by the high maximum lift coefficient for the power-on 
condition. Both of these phenomena illustrate the 
need for further study of the flow behind the propeller. 
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FIGURE 9.— Resistance-derivative coefTicient z„ of Doyle 0-2 airplane 

Uj angle of attack, degrees 
0 4 8 12 16 20 




Figure 10. — Resistance-derivative coefficient x u of Doyle 0-2 airplane 

The present investigation was made with n fixed 
stick, that is, a fixed elevator setting. This condition 
represents the case in which the stabilizer is set for 
the desired condition of flight and the pilot keeps the 
elevator in a constant position in spite of slight tend- 
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Figure 11.— Resistance-derivative coefficient z w of Doyle 0-2 airplane 
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Fig uit e 12.— Resistance-derivative coellicient x» of Doyle 0-2 airplane 
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Figure 13.— Rotary-derivative coefficient m w of Doyle 0-2 airplane 
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Figure 14.— Rotary-derivative coefficient m q of Doyle 0-2 airplane 
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Figure 15.— The period and damping coefficient for the power- 
off condition of Doyle 0-2 airplane 
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Figure 16. — Comparative curves of the measured period and 
damping coefficient for the power-off and power-on conditions 
of Doyle 0-2 airplane 
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cncies to fluctuate about the zero stick-force position. 
Another condition that might well be investigated is 
the case of free elevators as in that case the stability 
then attained determines whether or not the airplane 
can be flown hands-off . 

CONCLUSIONS 

The following conclusions are based on the results of 
the study of the dynamic longitudinal stability of the 
airplane used in this investigation : 

1. The theory of longitudinal stabilty based on the 
assumption of small oscillations is satisfactory for 
practical studies of the longitudinal stability of air- 
planes. 

2. The calculation of stability for the power-off 
condition from basic data is practicable and gives 



sufficient information to be of value in regard to the 
stability of new designs, though the results should not 
be interpreted too literally if they show the stability 
to be close to neutral. 

3. There is not enough knowledge at present con- 
cerning the effect of the slipstream on the various 
factors upon which the stability of the airplane de- 
pends to warrant calculations for the power-on con- 
dition. 



Langley Memorial Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Langley Field, Va., June 1932. 



APPENDIX 



DEVELOPMENT OF NONDIMENSIONAL STABILITY 
EQUATION 

The dimensional form of the equation of the longi- 
tudinal biquadratic, as given in the various works on 
stability, can be expressed as a determinant. 

m\ — X u -X w — XXf + W cos Bo+Xmwo 
— Z u m\ — Z w — \Z q +W sin d 0 — \mu 0 =0(1) 
-M u -M w + X 2 5-XM, 

where the subscript zero refers to the initial values of 
the variables and the derivatives are in terms of total 
force and moment rather than force and moment per 
unit mass. The nondimensional form was first sug- 
gested by the consideration 

X=-C D p%SV 2 

where the symbols have their usual significance. 

X u = dX/dV=-C DP SV 
XJpSV=-C D = x u 

where x u is the nondimensional form of the derivative 
X u . 

By similar reasoning, the coefficients of the other 
derivatives can be obtained and will be found to be 



Zu 

X q 
Za 



= ZJpSV 
= X W / P SV 
= Z W / P SV 
= XJ P SVl 
= Z Q / P SVl 
m u = M u /pSVlri 
m w = M u /pSVlrj 

m g = MJpSW v I 

where rj = B/ml 2 and I = a characteristic length. 

The choice of the length I is arbitrary. The distance 
from the center of gravity to the rudder post is con- 
sidered suitable for a study of longitudinal stability. 

Substituting for the dimensional derivatives in equa- 
tion (1) and replacing 

Wcos 6 0 by V 2 C L pSV 2 

Wsin 6 0 by y 2 C L pSV 2 tan 6 0 

the following is obtained 



and 



In any steady state, using the wind axes, w 0 will be 
zero and u 0 equals V. Then introducing the param- 
eter 

n = m/pSl 

the determinant can be written 
h — x u - x w ){C L \i - \x q 

— z u \ x — z w ) 2 C L p. tan e o — fx\ x — z a \i = 0 (2) 

— m u — m w X 2 ! — m Q \i 

In practice it has been found that x Q and z q are terms 
of the second order with respect to the other deriva- 
tives and consequently are hereafter neglected. 

The biquadratic obtained from the nondimensional 
determinant, which is of similar form to that obtained 
from the dimensional form, is 

\\ + B\\ + C\\ + D\ l + E=0 

where 

B= —m q — x u — z w 

(7= z w m q + z w x u + m Q x u — z u x w — p.m w 

D = y 2 jnm u C L + nm w x u + y 2 C L /xm^tan B 0 4- m q (z u x w 

) — X w \Ml u 

E=y 2 yiC L (x^m^-^mjtan d 0 + yip.C T Xm w z u -m u z w ) 

Stability is assured if Routh's discriminant R = BCD 
— D 2 — B 2 E and each of the coefficients are positive. 
For a more detailed analysis of the character of the 
motion following a disturbance, the complete solution 
of the equation is found. This analysis can be made 
by an approximate factorization of the biquadratic 
into two quadratics 

\ 2 1 + B\ 1 + C=0 

X 2 i + (^~^)xi + §=0 

of which the first represents a short-period, heavily 
damped oscillation, and the second a long-period, 
lightly damped oscillation. 1 It is with the lightly 
damped oscillation that instability is most likely to oc- 
cur and therefore it is usually only necessary to inves- 
tigate this phase. The solution for X x from the second 
quadratic is 



l/DC-BEV E 
\V 2C 2 J ~C 



BE\ 
2\C C 2 ) ± ' 

The motion is periodic only when the term represented 
by the radical is imaginary, and for this condition the 



-x wP SV 
m\—z w pSV 
—m w pSVlt] 



->a q pSVl+y 2 C L pSY 2j r\mw 0 



-\z q pSVl+y 2 C L pSV 2 tan 6 0 -\mu 0 



\ 2 ml 2 



-m q pSVl 2 r}\ 



m\—x u pSV 
-z u pSV 
—m u pSVlri 

Replace 

X by \ lP SV/m 
and cancel all factors common to every member of a 
line or column 



= 0 



1 The accuracy of this approximate factorization, which depends on the relative 
vaJues of the constants, was found satisfactory in the present case by reversing the 
operation with the appropriate numerical values of the constants substituted in 
the equations. 



Xi — x u —x w — XxXqpSl/m + C L /2 + \iW 0 jV 
-z u \i~z w - XxZgpSl/m + y 2 C L tan d 0 - \iU 0 /V 
— m u —m w — \ 2 x pSllm — \itn Q psllm 



= 0 
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period in nondimensional units is 

T 2tt 
1 IE / DC-BE \ 2 
\C A 2(7 2 ) 

The term — ^ — ^yP) is the constant part of the 

logarithmic decrement of damping and is called the 
damping coefficient f x . If the damping coefficient is 
negative, the motion is stable; that is, the oscillations 
are damped; and if the damping coefficient is positive, 
the motion is unstable. 

For comparison with the values obtained through 
the measurements in flight of the period and damping 
coefficient, it is necessary to convert the solution of Xi 
to the dimensional form. This conversion can be 
easily made by the introduction of the factor pVS/m. 
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Positive directions of axes and angles (forces and moments) are shown by arrows 



Axis 


Force 
(parallel 
to axis) 
symbol 


Moment about axis 


Angle 


Velocities 


Designation 


Sym- 
bol 


Designation 


Sym- 
bol 


Positive 
direction 


Designa- 
tion 


Sym- 
bol 


Linear 
(compo- 
nent along 
axis) 


Angular 


Longitudinal___ 
Lateral 


X 
Y 
Z 


X 
Y 
Z 


rolling 


L 
M 
N 


Y > Z 


roll 


0 

0 


u 

V 

w 


V 
<2 
r 


pitching 

yawing 


Z > X 

X > Y 


pitch 


Normal 


yaw 









Absolute coefficients of moment 

M 



Ct ~qbS 



Cm — 



gcS 



Angle of set of control surface (relative to neu- 
tral position), 5. (Indicate surface by proper 
subscript.) 



4. PROPELLER SYMBOLS 



Z>, Diameter. 

p, Geometric pitch. 

p/Dj Pitch ratio. 

V, Inflow velocity. 

V s , Slipstream velocity. 

T 

T, Thrust, absolute coefficient C T = 2 £p 



Q, Torque, absolute coefficient C Q = 



Q 



Pj Power, absolute coefficient C P = 

B lpV~ 5 

C s , Speed power coefficient = ^^* 
7) , Efficiency. 

n y Revolutions per second, r. p. s. 
3>, Effective helix angle = tan" 1 (j>~~) 



pn 2 D 5 

5. NUMERICAL RELATIONS 



1 hp = 76.04 kg/m/s = 550 lb./ft./sec. 
1 kg/m/s = 0.01315 hp 
1 mi./hr. = 0.44704 m/s 
1 m/s = 2.23693 mi./hr. 



1 lb. = 0.4535924277 kg. 
1 kg = 2.2046224 lb. 
1 mi. = 1609.35 m = 5280 ft. 
1 m = 3.2808333 ft. 



